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SUMMARY 

A  new  class  of  ionospheric  irregularities  in  the  auroral  oval  associated  with 
large  structured  plasma  flows  has  been  recently  isolated  with  radar  and  satellite 
in-situ  measurements .  These  density  irregularities  have  large  power  spectral  densities 
(pad)  at  short  scale  lengths  hundreds  to  tens  of  meters).  The  paper  characterizes 
the  density  and  velocity  spectra  in  such  regions  and  discusses  their  impact  on 
scintillation  observations  and  radar  performance,  The  structured  plasna  flows  nay 
occur  in  association  with  large  (MO  yA  m~z)  or  small  (M  y  A  m-z)  field  aligned 
currents.  The  velocity  spectra  have  fairly  shallow  power  spectral  indices  (*v-1.5)  in 
regions  of  large  field  aligned  currents  and  are  steep  (*v/-3)  in  regions  of  small  current 
flows.  The  density  spectra,  on  the  other  hand,  can  be  described  by  a  power  law  index 
^-2  in  both  the  large  and  small  field  aligned  current  regions.  The  temporal  structure 
of  scintillations  will  thus  be  dictated  not  only  by  the  scattering  strength  but  also  by 
the  large  flow  velocities  encountered  in  the  auroral  oval.  The  density  structures  with 
large  pad  at  short  scale  lengths  are  expected  to  introduce  considerable  ionospheric 
clutter  in  HF  radar  systems.  The  effect  of  E-region  conductivity  in  modifying 
small-scale  (<1  km)  P-region  structure  and  its  implication  for  communication  and  radar 
systems  is  discussed. 


INTRODUCTION 

In  recent  years  a  concerted  effort  has  been  made  to  study  high  latitude  F-region 
irregularities  and  their  relationship  to  bulk  plasma  processes.  Our  understanding 
improved  significantly  when  it  was  realized  that  km-scale  irregularity  structure  was 
associated  with  the  steep  edges  of  convecting  large  scale  plasma  density  enhancements 
both  in  the  polar  cap  and  auroral  oval  (see  Weber  et  al.,  1984;  1985;  Klobuchar  et  al., 
1986  and  references  therein).  Since  the  blobs  are  large  scale  features  on  the  order  of 
100  km,  they  have  long  lifetimes  and  can  be  tracked  by  incoherent  scatter  radars  and 
have  been  found  to  be  convecting  with  the  ambient  plasma  (Vickrey  et  al.,  1980).  The 
generalized  ExB  Instability  waa  found  to  be  a  viable  mechanism  for  the  small  scale 
structure  associated  with  the  blobs,  with  magnetic  field  aligned  currents,  neutral 
winds,  plasma  drifts  and  the  large  scale  blob-associated  density  gradients  all 
contributing  to  the  structure  formation  (see  Keskinen  and  Ossakow,  1983  for  a  review  of 
the  instability  processes).  If  radio  wave  propagation  techniques  are  used  as 
diagnostics  for  the  blob-related  irregularities,  then  one  would  see  scintillation 
structures  maximizing  on  the  edges  of  total  electron  content  (TEC)  enhancements 
(Basu  et  al • ,  1983) . 

More  recently,  another  class  of  F-region  irregularities  has  been  Isolated  which  is 
associated  with  large  shears  in  the  background  convective  plasma  flows.  The  unique 
feature  of  this  class  of  irregularities  is  that  it  is  not  associated  with  organized 
large  scale  gradients  of  plasma  density  (such  as  to  be  found  at  the  edges  of  blobs), 
but  rather  with  velocity  shears  with  shear  gradient  scale  lengths  of  1-10  km  (Basu 
et  al.,  1987a).  These  irregularities  are  also  found  to  have  large  pt>wer  spectral 
densities  (psd)  at  small  scales  (<1  km).  For  this  reason,  it  is  possible,  for 
Instance,  to  find  intense  scintillations  without  any  large  perturbations  in  TEC 
(Basu  et  al.,  1986). 

It  is  the  object  of  this  paper  to  present  simultaneously  observed  density  and 
velocity  spectra  obtained  in  the  vicinity  of  large  sheared  plasma  flows.  Lotova  (1981) 
has  discussed  the  effects  of  velocity  turbulence  on  the  modeling  of  scintillations  from 
in-situ  density  structure.  We  shall  compute  UHF  scintillations  and  HF  backscatter 
expected  from  such  density  spectral  fores  in  the  presence  of  considerable  velocity 
perturbations  and  discuss  the  effect  of  the  large  convective  flows  on  the  temporal 
structure  of  scintillations  observed  on  the  ground.  For  comparison,  we  shall  also 
compute  scintillations  using  observed  density  spectra  ij  high  conductivity  regions.  It 
is  found  that  in  such  regions,  for  example,  where  energetic  particle  precipitation  is 
observed  as  in  the  diffuse  aurora,  the  density  spectra  are  fairly  steep  with  little 
power  at  the  short  scales  (Basu  et  al.i  1984;  1987a).  The  reduced  power  at  such  scales 
(<100  m)  is  attributed  to  the  short  lifetimes  of  F-region  irregularities  In  the 
presence  of  underlying  E-region  conductivity  (Vickrey  and  Kelley,  1981;  Heelis  et  al., 
1985). 

OBSERVED  DENSITY  AND  VELOCITY  SPECTRAL  FORMS 

The  DE-2  satellite  (Hoffman  and  Schmerling,  1981)  la  an  excellent  source  for 
obtaining  high  resolution  F-rsgion  density  and  velocity  data  at  high  latitudes.  This 
data  base  has  recently  been  used  to  make  the  first  simultaneous  measurements  of  density 
(AN/N)  and  electric  field  fluctuation  (AE)  spectra  in  the  high  latitude  environment 
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(Basil  et  al.,  1987a).  It  is  to  be  noted  that  in  low  B  plasma  such  as  in  the 
ionospheric  F-region,  velocity  fluctuations  pre^electrostatic  and  equivalent  to 
electric  field  fluctuations  (i«e.#  v  »  E  x  B/B^).  Thus,  spectra  of  electric  field 
fluctuations  are  to  be  considered  equivalent  to  velocity  fluctuation  spectra. 

In  Figures  l(a»b)  and  2(a,b)  we  reproduce  3  samples  of  8s  each  density  and 
electric  field  data  and  their  respective  spectra.  The  symbols  S„,  Sg,  and  S.  on  the 
data  and  their  corresponding  spectra  represent  samples  taken  from  a  uniformly  high 
conductivity  region  (S„),  a  region  with  sharp  conductivity  gradients  (S^),  and  a  region 
of  low  conductivity  (S/).  At  the  time  the  DE-2  data  were  taken,  between  2149-2151  UT 
on  May  26,  1982,  the  satellite  was  traversing  the  dawn  sector  of  the  auroral  oval  (0700 
magnetic  local  time)  between  70-75°  invariant  latitude  in  the  southern  hemisphere  at  an 
altitude  of  approximately  350  km.  A  detailed  account  of  the  particle  precipitation 
character,  small  scale  field  aligned  Current  intensities  and  ion  and  electron 
temperatures  associated  with  these  samples  is  given  in  Basu  et  al.  (1987a). 

If  we  confine'  our  attention  to  the  density  and  electric  field  data  shown  in 
Figures  1(a)  and  2(a),  we  immediately  notice  the  lack  of  high  frequency  structure  in 
the  high  conductivity  region.  This  is  particularly  noticeable  in  the  electric  field 
data.  On  the  other  hand,  the  conductivity  gradient  region  (S«)  shows  a  very  large 
amount  of  structure  in  both  density  and  electric  field,  while  the  structure  in  the 
electric  field  is  somewhat  reduced  in  the  low  conductivity  region  (S,).  It  is 
interesting  to  note  that  regi  >n  Sp  is  associated  with  large  field-aligned  current 
intensities  of  tens  of  yA  m  in  Addition  to  the  large  velocity  shears,  whereas  S.  is 
associated  with  weaker  sheared  flows  and  smaller  field-aligned  currents. 

Before  we  discuss  the  (AN/N)  and  AE  spectra,  it  is  important  to  recognize  that  we 
had  to  make  the  assumption  that  measurements  of  spatial  structure  in  these  quantities 
made  on  the  satellite  are  Doppler  shifted  from  zero^to  the  frequency  f  -  kV  /2tt,  where 
V  is  the  satellite  velocity  (approximately  8  km  s  ).  In  this  paper  whenever  we  refer 
tS  wavenumber  k  or  wavelength  A  «  2ir/k,  it  is  based  upon  this  hypothesis.  Fortunately, 
there  is  convincing  evidence  for  this  assumption  of  "zero-frequency  turbulence"  (see 
review  by  Kintner  and  Seyler,  1985  and  references  therein).  Thus  we  feel  that  our 
assumptions  concerning  the  spatial  nature  of  the  density  and  electric  field  variations 
are  probably  valid  under  most  circumstances  in  the  ULF  frequency  regime  (Sugiura, 

1984),  except  possibly  when  propagating  waves  such  as  ion  cyclotron  waves  and  Alfven 
waves  are  encountered  over  limited  regions  (Temerin,  1978;  Gurnett  et  al.,  1984). 

The  spectra  of  (AN/N)  and  AE  in  the  three  different  regimes  are  shown  in 
Figures  1(b)  and  2(b)  respectively.  Both  (AN/N)  and  AE  spectra  are  fairly  steep  in  the 
high  conductivity  region  (S„),  while  both  are  fairly  shallow  in  the  conductivity 
gradient  region  (S_).  Moreover,  if  one  compares  the  power  spectral  density  (psd)  at 
10  Hz  (800  m)  of  the  two  density  spectra,  then  one  finds  approximatel y  15  dB  more  power 
in  the  S«  sample  than  in  S„.  In  the  AE  spectra  the  separation  in  psd  is  much  more 
dramatic  with  the  spectrum  showing  approximately  40  dB  more  power  at  5  Hz  (1.6  km). 
The  density  spectral  forms  and  psd  in  the  S.  and  Sg  regimes  are  almost  identical  with 
power  law  slopes  of  *v-2.  However,  the  slopes  of  tne  AE  spectra  in  the  S.  and  S_ 
regimes  are  quite  different  with  the  S.  spectral  slope  being  -3  comparedLto  thebS„ 
slope  of  shallower  than  -2.  In  all  3  situations,  the  important  point  to  note  is  that 
the  electric  field  (i.e.,  velocity)  spectra  are  either  equal  in  slope  or  steeper  than* 
the  corresponding  density  spectra.  Thus  following  the  arguments  of  Lotova  (1981),  it 
is  expected  that  we  should  be  able  to  model  scintillations  from  inysitu  density 
measurements.  However,  the  large  background  velocities  (100  mV  m""  is  equivalent  to 
2  km  8~  convection  velocities)  would  significantly  alter  the  temporal  structure  of 
scintillations. 


SCINTILLATION  ESTIMATES 

It  is  possible  to  estimate  scintillations  in  the  weak  scatter  formulation  (Rino, 
1979)  if  the  irregularity  amplitude,  density  spectral  forms,  background  densities  and 
irregularity  layer  thickness  are  known.  It  is  the  object  of  this  paper  to  present 
several  estimates  of  scintillations  that  may  be  observed  under  solar  maximum  or  minimum 
conditions. 

We  use  the  density  spectra  observed  in  the  velocity  shear  regions  to  estimate  the 
scintillations  to  be  expected  with  different  levels  of  background  ionization.  We  note 
from  the  top  two  panels  of  Figure  lb  that  the  psd  is  -18  dB  at  a  frequency  of  1  Hz. 
Using  the  satellite  velocity  of  8  km  s“  ,  this  frequency  is  equivalent  to  a  scale 
length  of  8  km.  We  use  the  measured  psd  and  the  observed  one-dimensional  density 
spectral  slope  of  -2  to  compute  intensity  scintillations  at  413  MHz  which  is  the  UHF 
frequency  transmitted  by  the  orbiting  HiLat  satellite.  We  utilized  the  observed 
irregularity  amplitude  of  approximately  20  percent  and  a  thickness  of  100  km  for  the 
irregularity  layer.  The  latter  seems  to  be  a  good  estimate  if  contours  of  maximum 
density  in  the  F-region  are  studied  using  latitude/altitude  scans  made  by  incoherent 
scatter  radars  at  auroral  latitudes  (Basu  et  al.,  1986).  The  weak-scatter  formulation 
of  Rino  (1979)  was  used  for  the  purpose.  The  equation  for  the  S,  index  (Briggs  and 
Parkin,  1963)  is  given  by: 

S42  ■  (reA)2  (L  sec6)  C$  Zv~>t  f(v)  J 


where  r  -  classical  radius  of  the  electron 
J  -  the  radio  wavelength 

L  -  the  thickness  of  the  irregularity  layer 
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0  -  the  ionospheric  zenith  angle 
Cs  -  strength  of  turbulence 
Z  -  Xzb  sec0/47r,  z  is  the  range 
v  -  related  to  l-Dspectral  index  given  by  2v-l 
J  -  geometrical  parameter  related  to  irregularity 
anisotropy 

f(v)  -  a  function  of  the  irregularity  spectral  index 
involving  T-functions. 


In  this  paper,  we  shall  confine  our  computations  to  the  S,  index  only  as  it  is  uniquely 
defined  by  the  measurement  of  Cg  and  V,  the  only  major  assumption  being  in  regard  to 
the  irregularity  anisotropy.  Tne  phase  scintillation  index,  on  the  other  hand,  depends 
on  the  detrend  filter  characteristics  and  the  effective  velocity  of  the  ray  path. 

Since  we  will  be  primarily  considering  auroral  observations ,  we  will  assume  sheet-like 
irregularities  aligned  with  the  local  L-shell  (Rino  et  al.,  1978).  The  best  estimates 
of  the  shape  is  given  by  8:4:1,  with  the  largest  anisotropy  being  along  the  magnetic 
field  (Fremouv  and  Secan,  1984),  The  effect  of  the  irregularity  anisotropy  on  the  S, 
index,  in  particular,  is  minimal  except  for  regions  of  exact  alignment.  4 

In  general,  the  greatest  degree  of  variability  is  observed  in  C<,  which  is  the 
strength  of  turbulence  (Fremouw  and  Secan,  1984).  The  in-situ  DE-2  data  thus  gives  us 
an  estimate  of  this  important  parameter  for  scintillation  computations.  As  noted 
earlier  the  measured  psd  of  -18  dB,  the  spectral  index  of  -2.Qand  irregularity 
amplitude  of  20  percent  translates  to  a  value  of  C_  -  4  x  101  mks  units  using 
equations  giv.en  by  Rino  (1979).  This  value  of  Cg  translates  to  S,  -  0.25  at  413  MHz 
for  nearQoverhead  locations  at  an  auroral  station  such  as  Tromso ,  Norwa y  (magnetic 
dip  -  78°).  If  field  alignment  is  considered,  then  the  S,  index  %  0.5.  Needless  to 
say,  at  such  times  the  S,  index  at  137  MHz  will  be  saturated  and  cannot  be  estimated 
using  this  weak-scatter  formulation. 

It  is  gratifying  to  note,  that  these  computations  agree  well  with  the  HiLat 
observations  at  413  MHz  made  at  Tromso  shown  in  Figure  3.  The  irregularities  during 
this  event  were  associated  with  velocity  shears  (Basu  et  al.,  1986).  Two  significant 
structures  were  observed  both  giving  rise  to  saturated  VHF  scintillations.  However, 
the  UHF  S,  index  is  slightly  less  than  0.3  for  one  and  0.4  for  the  other  which  is  in 
the  geometrical  enhancement  region. ^It  is  interesting  to  note  that  the  background 
density  at  this  time  was  3  x  101  ra-  as  compared  to  2  x  10  m”  used  in  the 

computations  based  on  DE-2  data.  This  indicates  that  large  irregularity  amplitudes 
coupled  with  relatively  shallow  slopes  can  cause  significant  UHF  scintillations  in  the 
auroral  oval.  Further,  when  the  background  density  is  larger  by  a  factor  of  3  or  4 
during  high  sunspot  conditions,  if  the  same  irregularity  amplitude  is  preserved,  it  is 
possible  to  have  a  three  to  four-fold  increase  in  the  scintillation  level.  Thus 
saturated  413  MHz  and  even  significant  1  GHz  scintillation  may  be  expected  in  the 
auroral  oval  in  conjunction  with  velocity  shears  with  particularly  intense  events  being 
observed  in  the  geometrical  enhancement  region. 


EFFECT  OF  UNDERLYING  E-REGION  CONDUCTIVITY  ON  SCINTILLATIONS 


The  bottom  panel  of  Figure  la  shows  density  data  (S„)  obtained  by  DE-2  in  the 
presence  of  diffuse  auroral  precipitation  which  creates  nigh  conductivity  in  the 
E-region.  The  field-line  integrated  Pedersen  conductivity  is  10  mhos  and  the  Hall 
conductivity  is  40  mhos  (Basu  et  al.,  1987a).  The  corresponding  spectrum  in  Figure  lb 
shows  the  spectral  slope  to  be  -2.5.  The  electric  field  spectrum  shown  in  Figure  2b  is 
even  steeper.  Figure  4  shows  spectral  data  in  the  auroral  oval  obtained  from  AE-D 
which  exhibits  a  one-dimensional  spectral  slope  of  -3  in  the  scalelength  range  smaller 
than  1  km  which  is  the  region  of  interest  for  UHF  scintillations  (Basu  et  al.,  1984). 

In  order  to  emphasize  the  effect  of  a  conducting  E-region,  we  choose  this 
experimentally  observed  slope  of  -3  in  association  with  the  large  irregularity 
amplitude  of  20  percent  which  was  observed  in  conjunction  with  the  velocity  shears 
discussed  earlier.  Using  exactly. the  same  background  density  and  anisotropy 
parameters,  we  obtain  Cg  ■  8  x  101  mks  units  and  S,  ■  .05  in  the  high  conductivity 
region  using  a  spectral  index  of  -3.  This  is  a  fairly  startling  result  which  shows 
that  scintillation  is  reduced  by  a  factor  of  5  for  the  same  irregularity  amplitude  and 
F-region  density  merely  by  going  into  a  region  whose  magnetic  field  lines  terminate  in 
a  conducting  E-region.  Vickrey  and  Kelley  (1982)  and  Heells  et  al.  (1985)  have  studied 
the  effect  of  E-region  conductivity  on  lifetimes  of  irregularities  and  we  point  out 
here  how  such  factors  can  have  a  large  impact  on  communication  systems.  The  upshot  is 
that  it  is  difficult  to  observe  significant  UHF  scintillations  in  regions  of  high 
E-region  conductivity.  This  is  an  important  factor  to  contend  with  in  the  auroral  oval’ 
and  polar  cap,  the  former  being  the  seat  of  the  energetic  auroral  precipitation  and  the 
latter  being  Illuminated  by  solar  radiation  continuously  in  the  summer  months.  Indeed 
intensity  scintillation  measurements  at  Thule,  Greenland  in  the  polar  cap  conducted 
over  half  a  solar  cycle  (1979  -  1984)  have  shown  the  enormous  decrease  of  scintillation 
occurrence  and  magnitude  during  the  local  summer  (Basu  et  al.,  1987b). 

EFFECT  OF  LARGE  VELOCITIES  ON  SCINTILLATION  STRUCTURE 


Figure  2  gives  examples  of  intense  small  sea le. ve] oci t y  structure  with  100  mV  m“ - - 

electric  field  giving  rise  to  a  velocity  of  2  km  s'*  in  the  high  latitude  environment. 

However,  at  other  times  fairly  steady  velocities  on  the  order  of  1  km  a”  nay  be 

observed  in  the  anti-sunward  direction  in  the  polar  cap  and  in  the  E-W  direction  in  Codd3 

auroral  latitudes  (Weber  et  al.,  1984;  Basu  et  al.,  1983).  Such  enhanced  convection  0/or 

velocities  have  a  noticeable  effect  on  the  temporal  structure  of  scintillations  as 
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manifested  by  measured  decorrelation  times,  intensity  spectral  shapes  and  the  intensity 
rate  of  scintillations.  We  present  intensity  scintillation  data  at  250  MHz  in  Figure  5 
from  two  consecutive  days  as  observed  at  Goose  Bay  within  the  auroral  oval 
(Basu  et  al.t  1985),  Significant  scintillations  of  the  same  level  were  observed  on 
both  days,  but  on  March  7  the  decorrelation  time  shown  in  dotted  lines  was  found  to  be 

2  s  while  that  on  March  8  was  0.4  s.  Estimates  of  the  velocity  could  be  obtained  from 
the  intensity  scintillation  spectra  at  the  2  times  denoted  by  arrows  on  Figure  5. 

These  two  spectra  are  shown  in  Figure  6.  The  spectrum  on  March  7  shows  a  Fresnel 
frequency  (fp)  of  0.05  Hz  while  it  is  0.3  Hz  on  March  8.  Using  the  relationship 

fp  =  u/(2AzP  where  u  is  the  drift  velocity  perpendicular  to  the  ray  path  (mostly. E-W 
fbr  this  geometry),  we  determine  this  drift  to  be  65  m  s“  on  March  7  and  390  m  s  on 
March  8.  Furthermore,  we  note  that  this  velocity  variation  gives  rise  to  widely 
different  intensity  rates.  For  the  data  sets  shown  in  Figure  5  the  intensity  rate  of 

3  ilB  s  occurs. in  0.1  percent  of  the  total  population  on  March  7;  while  the  rate 
becomes  7  dB  s  on  March  8  for  the  same  0.1  percent  level  as  shown  in  Figure  7 
(Basu  et  al.,  1987b).  Again  we  find  the  variation  of  a  magnet  os phe r i c  parameter, 
namely  its  convection  velocity,  having  a  significant  impact  on  propagation  systems 
operating  at  high  latitudes.  Phase  scintillations  were  not  discussed  in  this  paper, 
but  it  was  shown  by  Basu  et  al.  (1985)  that  enhanced  convection  yields  large  values  of 
phase  scintillation  when  the  source  is  a  geostationary  satellite.  Thus  convection  not 
only  affects  second  order  quantities,  such  as  decorrelation  time,  but  even  affects 
first  order  quantities,  such  as  phase  scintillation  magnitude. 


EFFECT  OF  UNDERLYING  E-REGION  CONDUCTIVITY  ON  F-REGION  RADAR  SYSTEM  OPERATING  AT  Hi¬ 


lt  was  mentioned  by  Basu  et  al.  (1987a)  that  one  of  the  unique  characteristics  of 
velocity-shear-associated  irregularities  was  the  existence  of  high  psd  of  (AN/N)  all 
the  way  up  to  approximately  10  m  scalelength.  In  Figure  8  we  show  the  psd  at  125  m  and 
46.5  m  observed  in  conjunction  with  the  three  categories  of  spectra  displayed  in 
Figures  1  and  2.  We  note  from  Figure  8  that  the  psd  associated  with  S^  is  almost  three 
orders  of  magnitude  less  than  that  observed  in  the  velocity  shear  region  indicated  by 
SG  and  S L.  The  filter  output  at  17  m,  while  not  shown,  indicates  a  similar  variation 
of  psd  between  the  two  categories,  that  in  the  high  conductivity  region  being 
10-**  Hz-*  whereas  it  is  10-®  Hz-*  in  shear  regions.  Now,  an  HF  radar  operating  at 
10  MHz  will  be  sensitive  to  a  scalelength  of  15  m.  Since  the  background  density 
variation  through  the  region  encompassing  S^,  S~  and  S,  in  Figure  8  is  insignificant, 
and  since  coherent  radar  backscatter  is  proportional  to  (AN;  (Woodman  and  Basu,  1978), 
we  can  expect  a  30  dB  decrease  in  radar  returns  from  the  F-region  at  10  MHz  when  the 
radar  beam  sweeps  into  the  high  conductivity  region  from  the  velocity  shear  region. 
Indeed  Baker  et  al.  (1986)  have  shown  that  HF  radar  backscatter  has  been  significant 
from  the  cusp  region  which  is  a  seat  of  large  velocity  shears  and  where  E-region 
conductivity  is  low  because  of  the  soft  particle  precipitation.  Thus,  we  find  that 
underlying  E-region  conductivity  controls,  to  a  great  extent,  the  evolution  of  F-region 
structure,  and  this  in  turn  has  significant  impact  on  radar  and  communication  systems 
operating  in  the  high  latitude  environment. 
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Intensity  scintillation  spectra  observed  on  March  7  and  March  8  at  tlaes 
denoted  bp  arrows  on  Figure  5.  Note  the  great  variation  In  Fresnel 
frequency  of  these  two  spectra. 
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